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We show that a generalization of the BCFW recursion relations gives a new and efficient method of 
computing correlation functions of the stress tensor or conserved currents in conformal field theories 
with an AdSd+i dual, for d > 4, in the limit where the bulk theory is approximated by tree-level 
Yang-Mills or gravity. In supersymmetric theories, additional correlators of operators that live in 
the same multiplet as a conserved current or stress tensor can be computed by these means. 



Introduction: The past few years have seen an in- 
tense study of the Britto-Cachazo-Feng- Witten (BCFW) 
recursion relations for gauge and gravity theories pQ. 
These recursion relations not only provide an extremely 
efficient route to S-matrix elements, they also do so with- 
out making any explicit reference to a local Lagrangian. 
Consequently, they are of interest both for computational 
reasons (such as the computation of S-matrix elements 
for physics at the Large Hadron Collider) and because 
they might help us shed new light on the formalism of 
quantum field theory. 

Hitherto, it has been believed that these techniques ap- 
ply only to quantum field theories in flat space. In this 
paper, we point out that the BCFW recursion relations 
can be generalized to compute sums of tree-level Wit- 
ten diagrams in quantum field theories in anti-de Sitter 
(AdS) space. By the celebrated AdS/CFT correspon- 
dence [2], this gives us new recursion relations for corre- 
lation functions in the dual conformal field theory (CFT). 

The calculation of Witten diagrams involving bulk 
gravity is notoriously difficult because of the infinite num- 
ber of interaction vertices that proliferate as the num- 
ber of external legs increase. As a result, even the four- 
graviton amplitude in AdSs (dual to the four-point cor- 
relator of the stress tensor in Af = 4 Super- Yang-Mills 
(SYM) at strong coupling) has never been computed di- 
rectly. The same problem exists in fiat space where 4D 
gravity has 2850 four-point vertices; the BCFW recursion 
relations reduce all computations down to the calcula- 
tion of the three-point function and greatly reduce this 
complexity. Our new recursion relations do the same for 
gravity in AdS, and we anticipate that they will simplify 
the computation of boundary correlators. 

The physical intuition underlying this study is as fol- 
lows. The BCFW recursion relations are predicated on 
the behaviour of Yang-Mills (YM) and gravity ampli- 
tudes when two of the external momenta are stretched 
off the infinity in a "complex direction." Although this 
is not strictly a high energy limit, it is nonetheless true 
that the amplitude is dominated by interactions between 
a soft background and a highly boosted particle at a sin- 
gle point. In this limit, we do not expect this highly 
boosted particle to see the curvature of the neighboring 
spacetime region. On the other hand, we do need to inte- 
grate over the different points where this interaction can 



occur. (This is similar to the intuition used in [3].) This 
process leads to the modified recursion relations that we 
present below. A higher-point correlator is broken down 
into the integral of the product of two lower-point corre- 
lators. Just as in fiat space, we can continue this process 
till we are left only with three-point functions. 

If we set out to compute a vacuum-correlator in the 
boundary theory, with all normalizablc modes switched 
off in the bulk, the recursion relations lead us to corre- 
lators computed in the presence of specific states; in the 
bulk this corresponds to turning on some normalizable 
modes. We will call these generalized correlators, transi- 
tion amplitudes. 

A further extension of our recursion relations allows us 
to compute transition amplitudes in supersymmetric the- 
ories, including Af = 4 SYM and the theory on multiple 
M5 branes in the supergravity limit. Perturbative com- 
putations in supersymmetric theories are often tedious; 
the recursion relations that we present ameliorate this by 
using a generalization of Nair's on-shell superspace [I]. 

Review of Perturbation Theory in AdS: We will 
work in Poincare coordinates where the metric is 

ds 2 = g^ u dx^dx v = z~ 2 {dz 2 + ijijdx' 'dx 3 j . (1) 

Poincare invariance in d dimensions makes it convenient 
to Fourier transform functions of x l and we will call the 
conjugate variables — ki — "momenta." 

For a non-interacting massless scalar field in AdS 
with spacelikc momentum, k, we have the unique solu- 
tion, 4> = foe ik - x z v K v {\k\z), where |fe| = |fc 2 | 1/2 and 
v = d/2, while for timelike momentum, we have the 
non-normalizable solution </> = foe %k ' x ' z v Y v i\k\z) , and 
the normalizable solution <f> — foe x z v J„(\k\z). We 
will use a uniform notation to write these solutions as 
cj> = foe lk ' x E v (\k\z), where z~ v E v is one of K v , J„, Y v . 

The free Yang-Mills equations in AdS are solved by 

Af(x, z) = e a E Vl (\k\z)e lk - x ; Al = 0, k • e a = 0, (2) 

where refers to the z-direction, v\ = v — 1, and the 
color index, a, is not italicized. We can move away from 
this gauge through A a (x, z) — > A a (x, z)+df Jl (/) a , where (j) a 
is any scalar field. Similarly, freely propagating gravity 
waves comprise transverse traceless tensors: 

hij = e tJ z~ 2 E„{\k\ Z y k - x :h 0fl = 0, he ij = 0, e\ = 0. 

(3) 
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We will refer to ef and e,j as "polarization vectors." 

The propagator for scalars, and for gauge-bosons and 
gravitons in this axial gauge is (J k = J ~ %d kdp 
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e ik < x - x 'h v - 2 J u (pz)J u (pz'){z'y- 
(k 2 + p 2 — ie) 

2-TijTki 
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where %j = rjij + kikj/p 2 and we have suppressed the 
trivial color dependence in G YM [5]. What will be im- 
portant for us is that, in each case, at p 2 = — fc 2 , the 
numerator of the integrand breaks up into a sum of a 
product of normalizable modes. 

Transition Amplitudes in AdS: Consider CFT 
operators 0(fe 31 ), . . . 0(k 3n3 ) and states s,s' that are 
dual, respectively, to linear combinations of normalizable 
modes with momenta fen, . . . ki ni and fc2i, ■ ■ • k 2ri2 hi 
the bulk. We will examine the transition amplitude 

T(k lm )(2ir) d 6 d (J2 k lm ) = (s\0(k 31 ) . . . 0(k 3n3 )\s'). 

Im 

Physically, we may think of |s');( s l as specifying data 
along the past and future horizons of the Poincare patch; 
we are then asking for the probability that the operators 
0(k 3m ) will induce a transition between these states. 

We will work at tree-level in bulk perturbation the- 
ory. To compute transition amplitudes, we draw bulk- 
bulk diagrams as usual. Then, we contract the legs with 
momenta in the set k 3m with bulk to boundary propaga- 
tors (non-normalizable modes), and the other legs, which 
carry momenta in the set k\ m or k 2m , with normalizable 
modes. The reader may prefer to think only in terms of 
this perturbative prescription and should consult [6] for 
further discussion. A vacuum correlator is just a special 
case of a transition amplitude, where all normalizable 
modes are switched off. 

The structure of perturbation theory tells us that 
transition amplitudes are produced by the action of a 
multi-linear operator on a set of (normalizable or non- 
normalizable) solutions to the equations of motion. For 
example, in Yang-Mills, with A^ (x, z) drawn from |2| 



T = G(A*\(x,z),...A%(x,z)). 
These operators obey Ward identities: 



(6) 



G(S7^ (x, z) 7 A% (x, z), . . . A^ (x, z)) = 0, (7) 
for any </> ai (x, z). An analogous identity holds for gravity. 



Below, we will consider transition amplitudes, 
T(fc m ,e am ), that depend on a set of discrete momenta, 
but also on polarization vectors for gauge-bosons and 
gravitons. The reader should note that some of the 
km, c am may correspond to normalizable modes, and oth- 
ers to non-normalizable modes; this will be left implicit. 

BCFW for Scalars: We start with a massless scalar 
cj) 3 theory to explain the main idea in a simple setting. A 
four-point transition amplitude involves three terms. 



T(fci,fca,fca,fc 4 ) 



E v {\k 1 \z 1 )E v {\kz\z 1 )z'iJ v (pz l ) 



x z v 2 J v {pz2)E v (\k z \z- 1 )E l ,{\k 4 \z 2 ) 



(fei +k 2 ) 2 +p 2 
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where the . . . are the t and u channel terms. Now, 
consider the extension kx — > kx + qw, fc 4 — > fc 4 — qw, 
which depends on the parameter w and where q 2 = 
q-fci = q ■ fc 4 = 0. This generically requires q to 
have complex components. Under this extension, the 
integrand of ^ is a rational function of w (although, 
of course, the integral itself is not) and has a pole at 
(2q ■ k 2 )w = —(p 2 + (fei + k 2 ) 2 ) with a residue 

[-lE^k^E^z^Mpz^- 1 ] 
x [-iz^J 1/ (pz 2 )E,{\k 3 \z 2 )E u (\k 4 \z 2 )z 2 d - 1 ] . 

Each bracketed terms is the integrand for a 3-point func- 
tion! There is also a pole at w — oo in ^ because the 
integrand of the diagram with a contact interaction be- 
tween kx and fc 4 goes to a constant at large w. 

It is easy to see that the same structure persists for n- 
point amplitudes. Poles in the integrand of a transition 
amplitude occur when the denominator of a propagator 
vanishes; the residue is the product of the integrands of 
two lower-point amplitudes and a simple factor from the 
propagator. The inclusion of the residue from w — oo 
permits us to completely reconstruct the integrand. So, 



T(fe a ,...fe T 



B 



E 

{7r} ,rr, 



-iT 2 



2{p 2 + K 2 ) 



dp 2 



(10) 



T 2 = T(ki(p), . . . Or(-C ■ • • Mp)). 



The sum is over all ways of partitioning the momenta 



' kn m } ; { k T 



1+1' • 



}, with 

in the other. Also, K = k± + 
K 2 +p 2 )/{2K ■ q);ki(p) = fex + 



into two sets {k±,k 
fei in one and fc, 
J2™ k^iwip) = 
qw(p);k n (p) = k n - qw(p);k' m = K - qw(p). The 
"boundary term," B, is the contribution from the pole 
at w — oo, comprising the sum of all diagrams where fei 
and k n meet at a point. 

Note that if we set out to compute a vacuum correlator, 
all the E in J8| are non-normalizable. Nevertheless, as 
we see from dOl, the mode corresponding to k' m in (10) 
will always be normalizable. This is implicit in (10). 
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Yang-Mills: Following [7], we expand the gauge field 
as A a = A a + a*. In background field gauge, the 
quadratic Lagrangian for is 

2C = D^alD> 1 a v > a + (2F^< a .f abc + R^S hc ) (11) 

where is covariant in spacetime and with respect to 
the background field; F is the background field strength; 
/ gives the structure constants and R is the Ricci tensor. 

We examine the large w behaviour of the two-point 
function for a®, BCFW extended as above. With q- 
lightcone gauge for the background field, q ■ A & = 0, all 
O (w) interactions come from the diagrams where there is 
a single interaction of the fluctuating field with the back- 
ground. This is because every propagator comes with a 
factor of (With the use of (|7|), we can show that 

the kikj terms in Q do not spoil this power counting.) 

From the effective action ( [IT] ), we see that the domi- 
nant contribution to the transition amplitude comes from 



:v7 D 



-2F 11 



",a b 
7 lM 



abc 



01A 

w 



d xdz 

■yd+l ' 



(12) 



where a!,a„ belong to ([2]). Below, we will suppress the 
color- factors, which are unimportant for our purposes. 
We choose the polarization for ai by e± = q, and define 



t by an 



- 1 (9 M - i M ) , where (f> = e^ kl+qu > x E Vl (z). 



By the Ward identity now, instead of ai^, we can use 
w -1 ^ in (12). As a result, the terms in the integrand of 



( 12 ) die off at large w if (a) e n does not grow at large w 
(which requires e n ■ q = 0) and (b) k± ■ e n = 0. In d = 4 
this forces us to take e n = q also. For d > 4, we can 
choose an e„ 7^ q that is orthogonal to fei, k n , q. 

With this choice of ei = q and these constraints on 
e n , we can reconstruct the integrand, up to terms that 
integrate to zero, using its poles at finite w. Repeating 
the argument above, we get the recursion relation 

-1T 2 
2(p 2 + K 2 ^ 



T(k 1 ,e 1 ,...k n ,e n ) = 

{7r} ,m,e' 



-dp 2 , 



T 2 = T(fci(p), e x , ...k' m , e' m )T(-k' m , e' m , . . . k n (p),e n ). 

(13) 

This has no boundary term and the sum now also runs 
over all normalized polarization vectors for k' m . 

These recursion relations are shown schematically in 
Figjl] Starting with a four-point vacuum correlator, we 
get the integral of the product of two three-point transi- 
tion amplitudes each of which has one normalizable mode 
(shown by the dotted line). 

It would be nice to find a way to relax our conditions 
on the polarization vectors. (In flat space, with ei = q, 
e„ can be arbitrary.) Nevertheless, by combining differ- 
ent BCFW extensions, we can compute the amplitude for 
arbitrary polarizations although we postpone the combi- 
natoric details to a forthcoming paper. 




FIG. 1. Recursion Relations 



Gravity: Expanding gravity fluctuations about the 
background metric, Q^ v = + /i M „, we find the gauge 
fixed quadratic Lagrangian [8 



S = 



-1 
64ttG 



d xdz 



p V + zhrR^hr} , (14) 



where W v = W - 



\g^h a P g a p, and all covariant deriva- 
tives are with respect to the background metric. 

It is now easy to extend the Ward identity argument 
above to gravity. We state the result. If we take the 
polarization vector for fei to be ejj = qiqj, then the inte- 
grand of amplitudes dies off at large w if we also take the 
polarization for k n to be either (a) = ^(iUj), where 



v ■ q = or (b) e-' 



V 



fei = 0. 



With this constraint, the gravity recursion relations 



v;,v^, where v m ■ q 



are the same as ( 13 1 with the obvious substitution of the 
gauge polarization vectors, e m , with gravity polarization 
vectors. We will describe the polarization-combinations 
that are accessible through combinations of BCFW ex- 
tensions in a forthcoming paper. 

Supersymmetric Theories: We can generalize 
these relations to supersymmetric theories as in flat space 
OHO]. AdSd+i supergroups do not exist for d > 6 and 
we are interested in the cases d = 4, 5, 6 [TT]. The reader 
may be more comfortable thinking about superconformal 
algebras in flat space although we will only use the super- 
Poincare subgroup of these algebras. A fact that we will 



use below is that ( 13 ) holds for Yang-Mills and gravity 



coupled to matter with the modification that the sum 
over polarizations must be expanded to run over these 
particles as well. 

We explain the case for d — 4 in detail and indicate 
results for d — 5,6. In d = 4, the superconformal group 
is 5/7(2, 2\Af). For J\f = 4, we have the 16 supercharges 
Qa,Qai an d their conformal partners. (We follow the 
conventions of [12 , so / is an R-symmetry index and a, a 
are spacetime spinor indices; see also [13].) We write 
the momenta we wish to extend, fei and k n , as linear 
combinations of two null vectors, A1A1 and A2A2, using 



(the a m i are some coeffi- 



cients) and take q a a — ^i a ^2a- Next, we assemble the 
vector of 27V-supercharges: = {(Q 1 , A2) , [Qi, Ai]}. 

(A runs over 1 . . . 2Af.) For Af = 4, defining T = 

r ^ixh1i^h r f li2 ' r f 1 ^i we nn d tna t; all states in the stress- 
tensor multiplet can be written as (with m = 1 or m = n) 



T m (v) = U + {n)T__{k m )U + {-n); U + ( V ) 



,Q+va 



(15) 
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The expansion of these operators in the 8 Grassmann 
parameters r\A contains all the original operators. With 
Af = 2, a similar expression exists for operators in the 
same multiplet as a conserved current. 

We pause to note that only half-Bogomol'nyi-Prasad- 
Sommerfeld (BPS) multiplets of the superconformal alge- 
bra can be represented using a form like (15 1. Although 



cial case, could be calculated by the recursion relations 



this form is not available for all half-BPS representations, 
it exists for all such representations in d = 4, 5, 6 that 
contain the stress tensor or a conserved current. 

Now, consider a n-point correlator that involves two 



operators from ( 15 1 with the same Grassmann parameter 
and n — 2 other operators, which we denote below by the 
composite operator Oq- The fact that this correlator is 
invariant under supersymmetry transformations implies 



(T^T^Oc) = (T__(fex)T__(fe n )0^ 



C/i 



(16) 



where 0' c = U+(—r])OcU+(ri). The right hand side can 
be computed by BCFW recursion as explained above. 

So, supersymmetry allows us to compute a "diagonal" 
subset of correlators i.e correlators of operators in the 
stress-tensor multiplet where at least two Grassmann pa- 
rameters are the same. This suffices to determine the full 
set of 4-point correlators in Af = 4 SYM, which can be re- 
duced to one independent function [2]. But, in general, 
we would like to compute correlators where all Grass- 
mann parameters are arbitrary. This is possible with 
flat space amplitudes; the difficulty here is that we have 
stricter constraints on the polarization-combinations that 
behave well under BCFW extension. 

In d = 6, the supercharges live in a 6 dimensional 
chiral-spinor representation (with eigenvalues ±1/2 un- 
der Lorentz transforms in the (2i — 1, 2i) plane) and in 
an R-symmetry group Sp(2Af) where Af is 1 or 2. (See 
[15] for conventions.) Apart from the "diagonal" sub- 
set above, d = 6 allows for another calculable subset 
of correlators: We choose fci = (1, 0, 0, 0, 0, 0), fc„ = 
(a, 6, 0,0, 0,0), 9l = (0,0, 0,0,1, /),<?„ = (0,0,1,7,0,0), 
and form two arrays of 4Af supercharges each: Qf- + = 
{Q±i/2,±i/2,i/2J' an d Qn+ — { ( 3±i/2,i/2,±i/2}- Then, 
for Af = 2, with T__(fc m ) = Ty(fe TO )^g4 and U m {rj) = 
cxp [Q„ + ?7a], we can compute any correlator of the form: 

(C/i(7 7l )?7„(ry„)T__(fe 1 )T__(fe n )O c C/ n (-ry„)C/ 1 (- 77 i)). 

(17) 

This is somewhat better than what we can do in d = 4. 

In d = 5, the supercharges are spinors under 5*0(5) 
and the R-symmetry SU(2). This algebra has a half- 
BPS multiplet containing a conserved current and we can 
compute diagonal correlators of operators in this multi- 
plet. However, the stress-tensor lives in a quarter-BPS 
multiplet |16) . So, not all operators in this multiplet can 



be reached via the analogue of ( 15 ) and we can only com- 



pute diagonal correlators among those that can. 

Results: We showed that transition amplitudes de- 
fined by |5|, which include vacuum correlators as a spe- 



( 13 1 for bulk Yang-Mills and gravity. Successive applica- 



tion of ( 13 ) allows us to relate all transition amplitudes to 



the three-point amplitude that is fixed, up to a constant 
factor, by conformal invariance. 

These relations are also applicable to interacting bulk 
scalars with the addition of a boundary term shown in 
(10 1. Supersymmetry allows us to compute additional 



correlators where we can convert at least two operators 
to conserved currents or stress tensors with appropriate 



polarizations. This includes the "diagonal" subset in ( 16 ) 
and, for d = 6, also includes operators of the form (17). 



In a forthcoming paper, we will apply these techniques 
to the calculation of higher order correlators. It would 
be very interesting to understand the deeper physical sig- 
nificance of these results and also extend them beyond 
tree-level in the bulk. 
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